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摘要 : 在 高 能 电磁 波段 ， 天 文 探 测 器 的 角 分 辨 率 较 低 ， 面 对 大 量 的 伽 玛 射 
线 观 测 数据 时 ， 研 究 人 员 需 要 利用 数据 分 析 软 件 迅速 找 出 一 些 未 知 天 体 的 位 置 。 
利用 2008 年 发 射 运行 的 国际 费 米 伽 玛 射 线 空 间 望 远 镜 〈) 的 高 能 伽 玛 射线 数据 
(大 于 兆 电 子 辐射 ，) ， 对 伽 玛 射线 暴 〈) ， 进 行 了 详细 的 快速 精准 定位 研究 ， 
主要 研究 了 不 同时 间 和 能 量 选择 时 对 天 体 源 的 探测 置信 度 的 影响 。 研 究 结果 表 
明 ， 选 择 伽 玛 射 线 暴 触发 时 间 零 点 〈) 到 以 内 的 时 间 范 围 可 以 很 好 地 定位 伽 玛 
射线 暴 的 高 能 对 应 体位 置 。 实 验 得 到 的 高 能 伽 玛 射 线 TS 位 置 图 最 佳 位置 与 低能 
电磁 波段 的 后 随 观 测 位置 很 好 的 符合 ， 表 明 实 验算 法 可 以 有 效 地 计算 该 类 天 体 
伽 玛 射 线 加 现 源 的 真实 位 置 。 
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近 几 年 来 ， 天 体 的 高 能 伽 玛 射线 波段 巡天 成 为 天 文学 研究 的 重要 手段 ， 完 
其 原因 在 于 伽 玛 射线 可 以 反应 各 类 高 能 物理 过 程 ， 如 宇宙 线 粒 子 加 速 和 辐射 过 
程 、 宇 宙 瞳 物质 作用 过 程 等 上 1 。 研 究 天 体 源 的 伽 玛 射线 能 谱 特 征 ， 可 以 证 论 某 
些 高 能 物理 过 程 是 否 主导 ， 如 中 性 on 介子 衰变 可 以 形成 峰值 流量 在 的 鼓 包 能 
谱 特征 由。 一 方面 ， 伽 玛 射 线 巡 天 可 以 有 目标 地 观测 (Target of 
Observation) 已 在 其 他 波段 有 观测 的 天 体 源 ， 如 哈 勃 望远镜 (Hubble Space 
Telescope) 的 光学 /红外 巡天 目录 ， 斯 隆 数字 化 巡天 (Sloan Digital Sky 
Survey, SDSS) 目录 以 及 澳大利亚 帕克 斯 天 文 台 64m 射电 望远镜 (Parkes Radio 
Telescope) 巡天 目录 加 ， 观 测 的 数据 可 进一步 限制 目标 源 的 物理 、 化 学 特性 。 
另 一 方面 ， 高 能 伽 玛 射线 探测 器 会 被 一 些 强烈 的 天 体 伽 玛 驾 现 源 触 发 ， 例 如 ， 
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存在 时 标 通常 大 于 能 区 的 时 标 〈 通 常 不 超过 100 FY) ， 在 暴 触 发 后 的 数 千 秒 也 


可 能 被 观测 到 中 5 ， 因 此 探测 器 可 应 用 于 探究 该 类 天 体 伽 玛 射线 暂 现 源 的 辐射 过 
程 。 随 着 我 国 天 文学 观测 的 发 展 ， 特 别 是 上 暗物质 卫星 (DArk Matter Particle 
Explorer, , B25) 、 硬 X 射线 调制 望远镜 (Hard X-ray Modulation 
Telescope, , SARS) 的 升 空 671， 它 们 的 科学 目标 集中 在 搜寻 一 些 未 知 源 的 
多 波段 辐射 并 进而 解释 一 些 极 高 能 的 物理 过 程 ， 这 些 过 程 一 般 在 实验 室 里 很 难 
实现 。 因 此 ， 快 速 定 位 末 知 伽 玛 射线 源 是 实时 跟 进 未 知 源 的 有 效 手 段 。 

由 于 在 高 能 伽 玛 波段 探测 器 的 角 分 辨 棕 较 低 ， 研 究 掌 握 精 准 定位 伽 玛 射 线 
源 的 算法 ， 对 于 伽 玛 射线 巡天 观测 起 着 很 重要 的 作用 。 实 验 研究 采用 国际 成 熟 
的 费 米 卫星 携带 的 大 面积 望远镜 (Large Area Telescope, ) 的 伽 玛 射线 巡天 
观测 为 实验 对 象 ， 选 取 伽 玛 射线 暴 为 暂 现 源 样本 ， 采 用 (1) 不 同时 间 选 择 ， 
(2) 不 同 中 心 位 置 选择 ， 计 算 伽 玛 射 线 暂 现 源 的 TS 位 置 图 ， 最 终 得 到 最 优化 
的 定位 。 其 次 将 得 到 的 位 置 ( 含 误差 ) 与 其 他 波段 的 后 随 观测 (Follow-up 
Observations) 比照 ， 调 整 实验 设 定 参 数 ， 从 而 得 到 实验 结果 与 实现 研究 目标 。 
1 费 米 卫星 巡天 观测 及 实验 样本 选择 

费 米 伽 玛 射线 太空 望远镜 (原名 ，) 是 多 国 合作 研发 的 ， 包 括 美国 、 法 国 、 
意大利 、 上 日 本 、 瑞 典 等 国家 。 它 的 科学 目标 主要 针对 活动 星系 核 、 脉 冲 星 、 伽 
玛 射 线 暴 、 宇 宙 线 等 。 搭 载 的 望远镜 视野 达到 全 天 区 的 20%， 于 2008 年 8 月 4 
日 开始 科学 观测 ， 并 于 一 年 后 实时 公布 观测 数据 。 伽 玛 波段 的 观测 能 量 范 围 从 
到 大 于 ， 单 个 入 射 的 天 体 伽 玛 光 子 被 各 层 探测 装置 逐 级 记录 ， 研 究 组 系 
析 了 这 些 光 子 事例 并 进行 分 类 ， 如 最 常用 的 Source fe Clean 类 型 光子 。 官 方 统 
计 从 2008 年 8 月 4 日 到 2011 年 8 月 4 日 ， 共 有 超过 1.8X10" 个 事件 ， 最 终 有 
约 1.44X10' 个 通过 更 细致 的 伽 玛 光 子 选择 标准 并 发 布 。 图 1(a) 是 5 年 观测 大 
于 伽 玛 射线 分 布 图 。 其 中 ， 银 盘 区 域 是 最 亮 的 伽 玛 射线 探测 区 域 ， 该 区 域 之 外 
存在 一 些 亮 点 天 体 源 ， 例 如 大 麦哲伦 星云 O, PREZ O.’ 


统 地 分 


(a) (b) 
E 1. (a) Sar AT ea NE 中 心 较 亮 的 区 域 是 银 盘 区 域 ; (b) 的 高 
能 光子 O 计数 率 随 时 间 分 布 图 
Fig. 1(a) Galactic Gamma-ray distribution observed by the Fermi-LAT in 5-yr 
mission, the bright part of which is galactic plane; (b) Count rates of high- 
energy photons as a function of time since GBM trigger (100 MeV-200 GeV) 
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年 10 月 30 日 。 图 1(b) 是 〈 数 字 为 触发 暴 的 命名 方式 ， 数 字 前 6 位 为 年 月 日 ， 
后 三 位 区 分 当天 观测 到 的 不 同 伽 玛 射线 暴 ) 的 光子 计数 率 的 时 间 变 化 图 。 从 图 
中 可 以 看 出 ， 在 触发 (Trigger Time, ) 后 100 秒 以 内 ， 光 子 计 数 率 呈现 脉冲 
式 分 布 ， 而 之 前 及 秒 后 ， 光 子 计数 率 比 较 平 稳 。 实 验 样 本 的 每 个 伽 玛 暴 ， 选 择 
到 的 的 光子 事例 。 样 本 选择 标准 为 在 触发 后 内 有 高 能 光子 被 探测 到 ， 即 能 量 大 
于 的 光子 。 通 过 初步 筛选 发 现 , 16 个 伽 玛 暴 ， 或 者 只 有 的 低能 伽 玛 光子 探测 ， 
或 者 在 较 晚 期 〈 之 后 ) 才 有 大 于 光子 探测 ， 因 此 排除 了 这 些 伽 玛 暴 样本 ， 共 收 
集 到 106 个 高 能 伽 玛 探测 的 伽 玛 暴 。 
2 Fermi 科学 软件 介绍 和 本 实验 算法 设计 
实验 是 基于 Fermi 官方 科学 软件 的 最 大 似 然 拟 合 分 析 。 首 先 ， 计 算 伽 玛 暴 
高 能 伽 玛 探测 的 test-statisticvalue (,， 当 达到 为 置信 度 ) 随 不 同 起 止 时间 
和 不 同 中 心 位 置 选择 的 分 布 回 。 单 次 数据 处 理 流程 包括 : 1) 光子 事件 类 型 、 
中 心 位 置 及 半径 (12°) 、 观 测 时 间 、 观 测 能 段 得 选择 〈 工 具 ) ; (2) 飞船 位 
置 相关 选择 〈 工 具 ) ; (3) 弥散 伽 玛 背景 计算 〈 工 具 ) ; (4) 光子 曝光 计算 
(LH); (5) RAMAMS CLR); (6) 置信 图 计算 (工具) 。 和 算法 设 
计 就 两 个 条 件 做 了 阔 值 条 件 设 定 ， 即 : 
1) 时 间 间 隔 选 择 。 对 单个 伽 玛 射线 暴 ， 所 有 的 时 间 起 点 〈) ABR He 
玛 暴 的 触发 时 刻 〈《) ， 而 时 间 终 点 〈《) 依次 累加 〈《) ， 如 时 ， 
分 别 为 () ，() #20) 。 


2) 中 心 位 置 。 该 步骤 中 有 两 个 位 置 。 (a) 初始 位 置 (Po) ， 以 的 初 
步 结果 为 初始 位 置 (Poo (b) 以 单 次 数据 处 理 流程 得 到 位 置 图 
的 最 大 值 的 位 置 (P) 再 次 执行 最 大 似 然 拟 合 ,得 到 本 次 数据 处 理 
流程 的 最 大 值 。 
基于 上 述 算 法 得 到 的 结果 画 出 特征 量 随时 间 选 择 的 分 布 图 。 这 些 量 是 : 
(a) 定位 的 定位 误差 半径 () ; (b) 最 大 TS 位 置 与 其 它 观测 位 置 距离 (O) ; 
(c) 每 次 最 大 值 〈) 与 比值 ( 记 为 ); (d) 相应 时 间 内 ， 预 估 来 自 伽 玛 暴 的 
AFAA O) 与 观测 的 总 光子 数目 O 的 比值 〈 记 为 ) o ERAF, FAKA 
动 得 到 所 有 选择 下 的 最 大 值 , 记 为 ， 此 时 的 终点 时 刻 为 。 
其 次 ， 选 取 上 述 的 时 间 范 围 ， 即 选取 ， 并 以 此 时 的 位 置 〈《) 再 执行 一 次 数 
据 处 理 流程 ， 得 到 最 终 的 位 置 图 ， 并 给 出 此 时 最 大 的 值 () 及 其 对 应 位 置 () 。 
最 后 ， 将 其 它 波 段 后 随 观 测 得 到 的 位 置 与 位 置 图 的 位 置 《) 及 其 置信 误差 
等 高 线 对 比 ， 测 试 件 玛 射线 源 位 置 的 准确 性 。 
3 实验 结果 及 讨论 
图 2 是 两 个 徊 玛 暴 的 、、、 等 由 设计 算法 得 到 的 特征 量 随 着 选择 的 截止 时 间 
O 的 分 布 关 系 。 蓝 色 竖 直 点 线 是 所 有 时 间 选 择 里 最 大 (BP) 时 间 位 置 《) 。 
从 图 中 可 以 看 出 ， 不 是 都 位 于 最 大 的 时 间 处 〈) ， 而 大 部 分 位 于 中 间 时 间 位 置 。 
这 可 能 有 以 下 3 个 原因 : (1) 在 之 后 ， 观 测 光 子 数 目 不 再 增加 ， (2) 在 之 后 ， 
所 在 的 飞船 位 置 受 到 其 它 因素 影响 ， 例 如 飞船 进入 南大 西洋 异常 区 域 (South 
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图 2. 两 个 伽 玛 暴 的 、、 随 不 同 的 变化 图 

Fig. 2Four observables as a function of selected end Time () 

Anomaly) 或 者 观测 中 心 轴 和 地 球 Limb K E, 使 得 河内 背景 过 强 ， 此 时 LAT 不 
能 正常 对 该 方向 计数 ; (3) 由 于 选择 的 光子 范围 为 12" 半径 的 圆 内 ， 在 之 后 的 
收集 的 光子 如 果 远 离 初 始 位 置 〈) ， 来 自 伽 玛 暴 的 光子 数 〈) 增加 较 小 ， 而 总 
光子 数 〈) 增加 较 快 ， 使 得 将 变 小 ， 从 而 会 使 得 值 变 小 。 另 一 方面 ， 当 定位 误 
差 半径 《) 和 位 置 距离 〈) 都 大 于 1° 时 ， 尾 应 大 于 +。 如 果 都 小 于 1 GN, m 

由 于 在 处 的 角 分 辨 率 (68% containment angle) 在 1° 左右 ， 因 此 不 要 求 必 须 
大 于 。 在 此 步骤 下 ， 得 到 合理 的 以 及 在 此 时 间 选 择 下 的 最 大 位 置 。 

图 3 是 在 选择 下 得 到 的 位 置 图 ， 该 图 给 出 了 两 个 暴 的 初始 位 置 〈《) ， 最 终 
RAZZE O ， 值 及 的 置信 误差 等 高 线 。 图 3 说 明 : (a) 这 两 个 暴 的 位 置 的 值 
都 大 于 25， 即 大 于 5 的 探测 置信 度 国 值 ; (b) GRB110731465 的 光学 波段 

(UVOT) 的 后 随 探测 位 置 《) 在 本 实验 定位 的 的 
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图 3. 两 个 伽 玛 暴 的 在 选择 下 得 到 的 位 置 图 。 三 角形 表示 初始 位 置 《) , UVOT. XRT% 
示 光 学 、X 射线 波段 的 后 随 观测 位 置 ， 十 字符 号 表示 本 实验 算法 所 定位 的 位 置 〈) ， 自 色 
等 高 线 是 的 置信 误差 距离 () 
Fig.3TS maps of two GRBs with the selected end time (). Triangle represents 
the initial position (), each of which is derived by UVOT or XRT labeled 


nearby. The cross is our resulting position (P3). The solid white line is 


contour line of localization at 95% confidence level. 

T 等 高 误差 线 以 内 ， 而 GRB120709883 的 射线 波段 的 后 随 探测 位 置 则 稍微 偏离 了 P3 
的 95% 等 高 线 内 部 ， 但 此 时 或 都 小 于 1 ， 考 虑 高 能 伽 玛 射 线 在 处 的 角 分 辨 率 C 
大 于 1*) ， 是 有 效 的 定位 。 

伽 玛 暴 在 keV 能 区 的 存在 时 标 ， 例 如 费 米 卫星 伽 玛 暴 的 〈Gamma-ray Burst 
Monitor, RX) ， 通 常 不 超过 100s， 比 较 了 在 时 间 内 和 在 选择 下 定位 高 能 
伯 玛 射线 的 位 置 的 精度 。 例 如 对 于 GRB110731465， 其 范围 为 ， 从 图 2 可 以 看 出 : 
时 间 下 的 为 8.6， 比 0- 时 的 280.4 小 很 多 。 因 此 ， 本 算法 能 够 精确 地 从 多 组 时 
间 段 里 找到 最 佳 时 间 段 ， 并 能 自动 校准 伽 玛 射线 中 心 位 置 以 确定 最 大 的 TS 位 置 。 


统计 的 探测 类 型 在 选择 下 最 大 的 分 布 情况 。 一 般 来 说 ， 当 时 ， 强 探测 置信 
R (Strong) ; 当时 ， 中 探测 置信 和 度 (Middle) ; 当时 ， 低 探测 置信 度 


(Low) ; 当时 ,无 明显 探测 (None) 。 结 果 显 示 的 样本 有 明显 的 伽 玛 射线 探测 , 
的 样本 有 强 的 伽 玛 射 线 探测 。 


表 1106 个 伽 玛 射线 暴 在 选择 下 的 分 布 情 
Table 1 Distribution of of 106 GRBs with the selected end time () 
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4 实验 总 结 

本 文 介绍 了 高 能 伽 玛 射 线 暂 现 源 的 位 置 探测 的 一 种 算法 ， 主 要 通过 探测 置 
信 度 〈) 随时 间 变 化 ， 考 察 了 不 同 蕉 止 时 间 对 探测 置信 度 的 影响 。 研 究 结果 表 
明 ， 截 止 时 间 通 常 位 于 到 之 间 ， 这 通常 可 能 由 于 飞船 运行 状态 、 中 心 轴 位 置 变 
化 或 者 本 征 的 伽 玛 射 线 光 子 分 布 造成 的 。 结 果 表 明 该 实验 方案 可 得 到 的 伽 玛 暴 
样本 的 明显 探测 和 的 样本 的 强 伽 玛 探测 。 综 上 所 述 ， 实 验方 案 可 快速 地 对 观测 
数据 进行 伽 玛 射线 源 的 位 置 的 确定 ， 该 位 置 可 通过 等 网 络 发 布 信息 给 其 它 空 间 
或 地 面 望 远 镜 、 探 测 器 ， 进 行 有 效 的 跟踪 观测 ， 从 而 对 未 知 天 体 伽 玛 射 线 源 进 
行 多 波段 研究 。 

本 实验 算法 可 作为 实验 室 开放 性 的 学 科 实 验 ， 对 计算 机 性 能 要 求 适中 ， 适 
合 具 备 天 文学 科 或 高 能 物理 等 硕士 点 以 上 的 院 校 进行 教学 和 科学 研究 之 用 。 在 
各 个 天 文 机 构 的 服务 器 中 ， 教 师 和 学 生 可 将 本 算法 作为 一 个 代码 包 应 用 于 其 他 
RR UE AR LIS i e e P, 
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Experimental research on a fast and precision localization on the 
astrophysical gamma-ray sources 
Qingwen Tang’ 

(1 Department of Physics, Nanchang University, Nanchang 330031) 
Abstract: With poor angular resolution of astrophysical detectors in 
the high-energy gamma-ray bands, the researchers are eager to a fast 
localization on the unknown sources, on which many development 


software will help deeply. In this work, we perform the details on the 


fast and precision localizations with different selections of time 
and energy spans on a sort of gamma-ray sources, Gamma-ray Bursts (), 
employing the high-energy observational data Q from the 
international telescope of the Fermi Large Area Telescope. The result 
shows that a selection from trigger time () to less than seconds 
after Q can localize the position of gamma-ray source detection 

The derived significance map of each GRB implies that the precision 
localization is consistent with the position at the low-energy band 
by the follow-up observations, which prove the validity of our 
experimental design. 

Keywords: gamma-ray sources; high-energy  astrophysicalsatellite; 


detection significance 


